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Abstract. As presented by Nice et al. (2000), long-term 
timing of the eclipsing binary PSR B1744-24A shows that 
the orbital period of this system decreases with a time- 
scaleof only ~ 200 Myr. To explain the much faster orbital 
period decay than that predicted by only emission of the 
gravitational waves (~ 1000 Myr) we propose that the or- 
bital evolution of this system is also driven by magnetic 
braking . If magnetic braking is to explain the rapid de- 
cay of the orbit, then A characterizing the effectiveness of 
the dynamo action in the stellar convection zone in the 
magnetic stellar wind formula must be equal to 1. 
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1. Introduction 

The formation and evolution of a very low-mass binary 
system consisting of a neutron star (millisecond pulsar) 
and a 0.1 M Q companion is not yet well understood. One 
of the main observational features of these systems is that 
they are eclipsing. One interesting eclipsing binary sys- 
tem, PSR B1744-24A, is in the globular cluster Terzan 5. 
The neutron star is an 11.56 ms pulsar in a 1.8 h orbit 
with a low-mass companion (Lyne et al. 1990). The du- 
ration of the eclipse in this system is very variable and is 
never less than one-third of the orbital period. Lyne et 
al.(1990) estimated that the energy flux from the pulsar 
at the companion surface (isotropic radiation is assumed) 
is < 2xl0 10 erg s _1 cm~ 2 which is near the critical value 
when the irradiation may influence the structure of a low- 
mass star (Podsiadlowski 1991). However, this value is at 
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least a factor of 15 less than for another eclipsing millisec- 
ond binary pulsar system PSR B1957+20 (Fruchter et al. 
1990). This small irradiative flux value and great varia- 
tions in duration of the eclipse lead Lyne et al. (1990) 
to the conclusion that in this system Roche lobe overflow 
may occur. 

Recent long-term timing observations of the eclipsing 
binary pulsar PSR B1744-24A at the VLA and Green 
Bank, show that the orbital period of PSR B1744-24A 
has decreased, with a time-scale of only P or b / 'P or (,~200 
Myr. This is five times faster than if the orbital period 
decay is driven only by the gravitational wave radiation 
(Nice et al. 2000). 

The key question of interest is: how did this system 
evolve? 

In this short note we show that the new observational 
data for the orbital period decay of PSR B1744-24A may 
give us a unique opportunity to test the magnetic braking 
mechanism if the observed orbital period decay is secular. 

2. General picture of low mass binary system evo- 
lution 

More than twenty years ago it was realized that to produce 
short orbital period low-mass binary systems (cataclysmic 
variables and low-mass X-ray binaries - LMXB) it is nec- 
essary to include two mechanisms for orbital angular mo- 
mentum losses: gravitational wave radiation and/or mag- 
netic braking. Low-mass unevolved stars (main-sequence 
stars) will fill their Roche-lobe only for orbital periods less 
than 12 hours and the mass transfer due to gravitational 
wave radiation mechanism never exceeds ~ 10~ 10 M Q /yr. 
This produces an X-ray luminosity of ~ 10 36 erg/s. Obser- 
vations, however, show that many LMXBs in this orbital 
period range have X-ray luminosities that are one or two 
orders of magnitude higher. 

To explain this discrepancy Verbunt & Zwaan (1981) 
introduced additional orbital angular momentum loss by 
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magnetic braking which was able to drive a much higher 
mass-transfer rate. In view of these much higher mass 
transfer rates the companion will be driven out of ther- 
mal equilibrium. However magnetic braking will not work 
forever. As soon as the companion mass has decreased 
to M ~ 0.3 M Q (fully convective star) an orbital period 
near the 3 hours, the rate of angular momentum loss by 
magnetic braking is expected to either vanish or to drop 
considerably (Spruit & Ritter 1983), causing the mass- 
loss rate to drop, and the secondary star subsequently 
relaxes to its thermal cqulibrium state. The stellar radius 
decreases and the star detaches from its Roche lobe. Dur- 
ing the detached phase, orbital angular momentum loss by 
gravitational wave radiation will continue, and when the 
orbital period has been reduced to 2 hours, the companion 
(of mass ~0.3 M Q ) again fills its Roche lobe. This is the 
so-called standard cataclysmic variable (CV) evolutionary 
scenario. 

In 1985 Tutukov et al. discussed what will happen with 
a low-mass binary if the star filling its Roche lobe is a 
slightly evolved star (so-called "turn-off main-sequence 
star"). In this scenario the binary evolves towards very 
short orbital periods. Due to the chemical composition 
gradient, the secondary star does not become fully con- 
vective when its mass has reduced to 0.3 M Q and mag- 
netic braking does not vanish. Between orbital periods of 
two to three hours, the mass transfer rate is very low, 
and apparently during this time the system may be in 
the "propeller" stage (Ergma & Sarna 1996). According 
to the standard CV evolutionary picture, the magnetic 
braking mechanism vanishes near an orbital period of 3 
hours. However in the latter scenario, it will also work for 
smaller values of the orbital periods . 

3. Model 

The evolutionary sequences for the secondary star in a 
low-mass binary were computed using a standard one- 
dimensional stellar evolution code based on a Henyey type 
code developed by Paczyhski (1970) and adapted to low- 
mass main-sequence stars. For more detail about the com- 
puter program see Muslimov & Sarna (1993) and Sarna & 
De Greve (1994, 1996). 

We compute the loss of orbital angular momentum due 
to gravitational radiation using the formula presented by 
Landau & Lifshitz (1971). 
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where Mi, M 2 , M and P oro are the primary mass, sec- 
ondary mass, total mass and orbital period of the system. 
We assume that the companion is being spun down by the 
magnetic braking and its spin and orbital rotation are syn- 
chronized at the cost of orbital angular momentum loss. 
We adopt the standard formula (sec Mestel 1968, Mestel 



& Spruit 1987) for the rate of orbital angular momentum 
losses due to magnetic braking of the companion. 
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where k 2 is the radius of gyration of the secondary 
star, k 2 =0.1. This formula contains a poorly known pa- 
rameter A (~ 0.7-1.8) characterizing the effectiveness of 
the dynamo action in the stellar convection zone (Verbunt 
& Zwaan 1981). To take account of the angular momen- 
tum loss that accompanies mass loss during the "propeller 
phase" , we use a formula based on that used to calcu- 
late angular momentum loss via a stellar wind (Paczyhski, 
1967, Sarna & De Greve, 1994) 



J 

J WIND 



hh 



M X M 2 
M 2 M ' 



M = /iM 2 , Mi = -M 2 (l-/i)(3) 



where j\ 



M acc 



and M 2 (<0) is the mass-loss rate 



from the secondary star; Mi(>0) is the accretion rate on 
to the neutron star, M(<0) is the rate of the total mass 
loss from the system; fx is the ratio of the mass ejected by 
the neutron star to that accreted by the neutron star and 
f 2 is defined as the effectiveness of angular momentum 
loss during mass transfer (Sarna & De Greve, 1994,1996). 

We consider the following model (similar scenarios 
have been discussed by Shaham & Tavani 1991, Kluzniak 
et al. 1992 and Ergma & Sarna 1996) for the PSR B1744- 
24A system. The progenitor of this system is a low-mass 
(<1-1.5M ) star + old neutron star. The secondary star 
fills its Roche lobe, as a "turn-off main-sequence star" . 
At the beginning, the mass transfer rate proceeds on a 
thermal time-scale and the neutron star will spin up to 
a millisecond period. After that, the mass transfer rate 
drops and the system enters the "propeller stage" and the 
pulsar spins down. 



Table 1 



Case 


A 


Port/10- 12 


T 






[ss- 1 ] 


[Myr] 


(A) 


0.7 


1.16 


174 


(A) 


0.9 


1.10 


184 


(B) 


0.7 


2.44 


83 


(B) 


1.0 


0.83 


244 


(B) 


1.2 


0.69 


292 


(B) 


1.5 


0.60 


337 


(C) 


0.7 


1.77 


114 


(C) 


0.9 


1.21 


167 



As an example we calculated several evolutionary se- 
quences: (I) M 2 = 1M Q , Mi=1.4 M , P;(RLOF)= 1 day; 
(II) M 2 = 1.5M , Ml=1.4 M q P;(RLOF)= 1.02 days 
(Pi(RLOF) is initial orbital period when the secondary 
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star fills its Roche lobe). In Fig. 1 (a, b, c) the mass- 
loss rate versus orbital period is presented for the follow- 
ing cases: (A) sequence (I), /i=/2=0, (B) sequence (I), 
/i=/2=l, (C) sequence (II), /i=/2=l- At first, the mass 
exchange rate is rather high (proceeds on the thermal 
time-scale of the secondary star) and during this phase, 
the old neutron star will spin-up. Latter M decreases and 
near an orbital period of ~ 2-4 hours, the mass accre- 
tion rate has its minimum value. Depending on the value 
of the surface magnetic field strength, the system may or 
may not be in the "propeller phase" . How does the value 
of A influence the orbital evolution of the binary system? 
For example for sequence (A) and for A= 1.8 , the final 
orbital period is larger than the initial and for A =1.0 a 
short period binary system with low-mass helium white 
dwarf and millisecond pulsar is formed (P or b ~ 13 hours). 
Only for A < 1 does the system evolve towards very short 
orbital periods (Fig.l). 
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Fig. 1. The dependence of the mass loss rate as a function 
of orbital period for three different evolutionary sequences (A) 
(upper panel), (B) (middle panel), (C) (lower panel). For de- 
tails see text. 
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Fig. 2. The evolution of the orbital period derivative with the 
orbital period. Other details the same as for Fig. 1 



In Fig. 2, a, b, c the variation of P or b versus orbital 
period is shown for the same cases as presented in Fig. 
1 a, b, c. As is clearly seen in this Figure, orbital period 
changes depend on the value of A. In Table 1 we present 
P or b and orbital period decay time-scale T=P or b / P or b val- 
ues for various A when the orbital period is equal to 1.8 
hours. The evolutionary time-scale obtained by Nice at 
al. (2000) is 200 Myr. From Table 1 we can see that a 
A close to 1 agrees well with the observed value for PSR 
1744-24A.The calculated mass of the secondary ~ 0.15 
Mq near a period 1.8 h also agrees well with the observed 
value (Lyne et al. 1990). 

4. Conclusion 

If the observed orbital period decay is secular and our 
suggestion is correct that the orbital evolution of PSR 
1744-24 A is driven by magnetic braking then this may be 
the first case when the orbital period decay by a magnetic 
braking mechanism has been measured. According to our 
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model calculations the best fit to the observed orbital pe- 
riod decay is obtained when A is near to one. Which means 
that the value proposed by Smith (1979) A=1.78 underes- 
timates, that and proposed by Skumanich (1972) A=0.73 
sligtly, overestimates the magnetic braking efficiency 
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